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Abstract 

X-ray luminosity function derived from observations of X-ray sources 
in galactic bulges can be explained by prinicipal evolutionary relations 
for mass accretion rate onto the compact object. The observed mean 
distribution of individual X-ray luminosities of galactic LMXB is sat- 
isfactorily described by a symmetric quasi-Lorentzian curve. The flux 
variance for bright sources is found to be proportional to the mean lu- 
minosity. Such a distribution does not change the slope of the power 
law luminosity function of the source population, which is expected 
from the dependence of the mass transfer rate on the mass of the 
Roche-lobe filling non-degenerate optical component. 



1 Introduction 



X-ray flux measurements from numeorus point-like sources in other galaxies 
have recently become possible by Chandra and XMM-Newton X-ray tele- 
scopes. Most of these sources are accreting compact stars in binaries. These 
observations have important implications for population studies of X-ray bi- 
naries in galaxies. Grimm et al. (2993) and Gilfanov (2004) carried out a de- 
tailed statistical analysis of the observed X-ray luminosity function of X-ray 
sources observed by Chandra and XMM. They constructed X-ray luminosity 
function of X-ray binaries in individual galaxies, derived the mean luminosity 
function and discussed its astrophysical implications. In our previous paper 
(Postnov 2003) we noted that the universal power-law shape of X-ray lumi- 
nosity function (XLF) of high-mass X-ray binaries obtained by these authors 
dN / dL x ~ L~ L6 in a wide range of X-ray liminosities 10 34 — 10 40 erg/ s can be 
readily explained by the universal properties of accretion from stellar wind 
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from the early-type component. Important points were the possibility to use 
the power-law dependence for mass-luminosity and mass-radius relations for 
early-type stars, as well as a natural assumption on the power-law initial 
stellar mass function (Salpeter law or Miller-Scalo law). We also assumed 
that the mass accretion rate onto the comact star is the principal parameter 
determining the X-ray luminosity: L x oc M a . 

New observational data were used by Gilfanov (2004) (see also Kim and 
Fabbiano 2004) to construct XLF of low-mass X-ray binaries (LMXB). It 
turned out that the mean XLF for these sources, too, can be approximated 
by a univeersal function, which is, however, not a single power-law, as in the 
case of HMXB. 

It is well known that X-ray binaries in general are variable. Moreover, many 
of them are transient sources, with X-ray luminosity changing by several 
times and even orders of magnitudes (e.g., X-ray novae) over different time 
scales. Typical X-ray observations do not last for more than several tens 
kiloseconds, which in most cases less than the characteristic time of signif- 
icant luminosity changes in individual sources. So the question arises as to 
what exactly the observed mean XLF reflects, the instantaneous luminosity 
distribution of a collection of sources or the luminosity function of individual 
sources. 

In the present paper we first show that the observed dispersion of individual 
X-ray luminosities of galactic LMXB as obtained from the analysis of the 
RXTE ASM data can be ona average described by a symmetric Lorentzian 
distribution with the variance proortional to the mean flux. Nest we prove 
that such a distribution of accretion rate M a onto the compact object does 
not change a power-law character of XLF which is theoretically calculated 
from the analysis of mass exchange rate M Q in close binaries. We arrive at 
the conclusion that the mean XLF of LMXB can be explained by evolution- 
ary features of the mass transfer process in low-mass close binary systems: 
the mass transfer rate due to binary's angular momentum removal by gravi- 
tational waves shapes the low-luminosity end of of the observed XLF oc L~ x 
(L x < 2 x 10 37 erg/s), while mass exchange rates driven by magnetic stellar 
wind from the late-type optical component can reproduce the high-luminosity 
end of the observed XLF of LMXB oc L~ 2 (2 x 10 37 <L,< 5 x 10 38 erg/s). 
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2 Power-law luminosity functions 



First of all, let us show on general grounds how power-law XLF in accreting 
binaries can be obtained. The basic assumption is the possibility to express 
the X-ray luminosity of a source through the mass exchange rate M Q which 
is a powe-law function of the optical star mass M Q (Postnov 2003): 

L x ~ M a ~ M ~ M: (1) 

Below we shall consider the luminosity function per logarithmic interval 
dN/d\nL x . Under assumption (1) in the stationary situation we have: 

dN dN dlnM _2«t=l 

- oc L x a (2) 



d\nL x d\nM dlnL 
Here the stationary mass function for optical components was utilized: 

(for example, (3 st = 2.35 for the Salpeter mass function). Consequently, a 
power-law shape for XLF is recovered: 

w r 



X 



d\nL x 

where the power-law index is 

r _ Pst - 1 

a 

and depends only on the slope of the mass function and relationship (1). 



2.1 Stationary mass function of optical components 

Two cases should be distinguished: The case of young massive close binary 
systems in which the mass of the optical component virtually did not change 
from the original value, and teh case of old low-mass close binaries, where 
the mass of the optical star significantly diminished in the course of mass 
transfer. In the former case the mass fucntion of the optical component does 
not alter significantly, so we can adopt its initial form: 

f(M) st = f (M) oc M-P f3 st = (3 (4) 
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In the latter case the stationary optical star mass distribution can be derived 
from the one-dimensional kinetic equation with stationary source f Q {M) 
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f(M)M oc/ (M) (5) 



DM 

In our case M < and we obtain 

• , s [m™* UM')dM' 

MM) = ^ (6) 

where M max is the upper mass limit for stars. For power-law mass fucntions 
its exact value is of minor importance. 

Substituting the power-law mass function f (M) ~ M _/3 into Eq. (0), we 
find 

f st (M) oc M-P +1 ~ a , (7) 

that is the power-law index in the stationary stellar mass distribution be- 
comes dependent on the mass exchange rate: 

(3 st = f3-l + a (8) 

For example, for HMXB accreting from steallr wind a ~ 2.5 (Postnov 2003) 
and for j3 = 2. 35. ..2. 5 we get T = 0.54. ..0.6 (the observed value is ~ 0.6). 



3 Effect of the luminosity function of indi- 
vidual sources on the luminosity function 
of the source population 

Before considering a more complicated case of low-mass close binary sys- 
tems, let us discuss the influence of individual source XLF on the XLF of 
the entire population. Let the accretion rate onto the compact object and 
hence X-ray luminosity L x of a source changes according to some distribution 
function F(x) (i.e., F(x) gives the probability to detect the X-ray luminosity 
in the range [F(x),F(x + Ax)]). Here the mean mass exchange rate M a is 
allowed to change or remain constant. Let Lq be the X-ray luminosity cor- 
responding to the mean accretion rate. We assume that this rate represents 



4 



Table 1: Low-mass X-ray binaries with short orbital periods. For each source 
given are the orbital period P(h) (in hours), the mean ASM flux F x (in ASM 
count rate), the rms variance of the flux o (in ASM count rate), minimal 
(F min ) and maximal (F max ) ASM fluxes over the entire time interval used in 
the analysis 



N% 


Source 


P(h) 


F x 


a 


F ■ 

± mm 


± max 


1 


V1333 Aql 


18.97 


2.31 


7.85 


-13.85 


49.68 


2 


V821 Ara 


14.86 


6.80 


13.51 


-13.92 


72.57 


3 


V1408 Aql 


9.33 


1.89 


1.23 


-17.07 


6.95 


4 


LU TrA 


9.14 


1.01 


1.06 


-4.59 


8.63 


5 


V691 CrA 


5.57 


1.18 


1.33 


-8.62 


16.73 


6 


V926 Sco 


4.65 


12.89 


3.38 


-9.70 


29.64 


7 


V2216 Oph 


4.20 


17.34 


2.51 


0.00 


35.78 


8 


GR Mus 


3.93 


1.93 


1.04 


-6.02 


23.36 


9 


V801 Ara 


3.80 


10.90 


4.88 


-5.63 


34.92 


10 


1705-4402 


1.31 


11.67 


8.85 


-3.89 


35.94 



the mass exchange rate between the components determined by the system's 
evolutionary state. It is easy to see that the observed XLF is the convolution 

f^ max (dNl dL )F(L x - L )dL 

dN/dL x = JLmin \ ' — —— (9) 

F{L X - L )dL 

We restrict ourselves by considering power-law XLFs, i.e. dN/dlnLo ~ Lq T . 
If the function F(x — y) can be factorized in the form 

F(x-y) = F 1 (y)F 2 (x/y-l), (10) 

then the power-law XLF of the source population dN/dlnx ~ x~ r pre- 
serves when Fi(y) is also a power-law function for arbitrary function F 2 (all 
necessary integration conditions are assumed). The proof readily follows 
from changing variables t = x/y — lin integrals and setting L m i n — > 0, 

Lmax * 

Let us use the RXTE ASM data to find the shape of the function F(L X — 
Lq) for galactic LMXBs. First let us consider bright LMXBs with known 
orbital periods shorther than 20 hours and mean fluxes above 1 ASM count/s 
(Table 1). In close binasries with orbital periods below 20 hours the mass 
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Figure 1: XLF for 10 LMXBs from Table 1 (with known orbital periods short- 
her than 20 hours) constructed from the RXTE ASM data. The luminosity 
of each source is given in units of the mean luminosity L . Distributions are 
normalized to unity J ^-dL x = 1. 
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Figure 2: The mean XLF for 10 LMXBs from Table 1 (the histogram). The 
solid lines shows the approximation by a Lorentzian curve. The confidence 
level is P(xj 7 > 8.3) « 0.96. 



of the optical star filling Roche lobe is < 2M (we assume the compact star 
mass 1.4 M and main-sequence optical star with solar chemical abundance). 
The mass transfer is driven by the orbit angular momentum removal due 
to magnetic stellar wind or gravitational wave emission (for orbital periods 
shorter than several hours, see below). Table 1 lists the mean flux of the 
source F x (averaged over several years of observations, typically 10 3 — 10 4 
individual points per source) and rms deviation a from the mean F x . Both 
F x and a are given in units of the ASM count rate (we remind that the ASM 
flux from Crab is 75 ASM counts/s). 

In Fig. 1. we plot indivudual normalized XLF of sources from Table 1 
calculated with the use of the ASM data and show the values of the corre- 
sponding orbital periods. Fig. 2 shows the mean XLF calculated for these 
systems (the histogram) and its approximation by a quasi-Lorentzian curve 
F{x — y) — ( x _ y y2 +by 2 (a, b = const), which obviously satisfies condition 
(|l()jl . The errors in each bin are calculated using the formula for variance 
(jj = \J~Y^j=\ ( x j ~ %)/ n i i n i is the number of points int he bin, x is the mean 
value). The formal confidence level of the fit according to the x 2 criterion 
for 17 degrees of freedom is P(xi7 > 8-3) ~ 0.96. 

As is seen from Table 1, almost all sources with known orbital periods are 
weak, so individual XLFs (especially left wings corresponding to weak fluxes) 
are determined with substantial errors. So as an additional test we used a set 
of bright LMXBs with arbitrary (or unknown) orbital periods (Table 2). XLF 
of an accreting X-ray source is determined by properties of non-stationary 
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Table 2: Bright LMXBs. The mean fluxes F x and the rms flux variances a 
are given in units of the ASM count rate. 



N% 




r X 


rr 

U 


1 
1 


byg A2 


o 1 .22 


Q QO 
0.62 


2 


OCO At 


oon one; 


11/1 fifi 

114. Oo 


Q 

o 


LtA y+y 


1 1 .22 


O A O 


4 


LrA y+i 


on o/i 


4.0 / 
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6 


GX 3+1 


22 34 


5 99 


7 


GX 17+2 


44.86 


5.33 


8 


GX 13+1 


22.45 


2.93 


9 


GX 340+0 


29.39 


4.34 


10 


GX 349+2 


51.40 


12.19 


11 


V926 Sco 


12.89 


3.38 


12 


V2216 Oph 


17.34 


2.51 


13 


V801 Ara 


10.90 


4.88 


14 


1705-4402 


11.67 


8.85 



accretion on the compact object. The external radius of the accretion disk 
is moderately depends on the system's orbital period (a oc P 2 / 3 ). The mean 
flux is determined by the mass of the compact object. In our case compact 
objects are neutron stars of about the same mass. Individual XLFs of bright 
LMXBs from Table 2 are shown in Fig. 3, the mean XLF for these sources is 
presented in Fig. 4 (the histogram). As for LMXBs with short orbital periods 
shown in Fig. 2, the mean XLF for bright LMXBs is well approximeted by 
the Lorentzian curve (confidence level P(xh > 8-8) ~ 0.94). 

Fig. 5 demonstrates the variance of the observed XLF for birght sources 
from Table 2 as a function of the mean flux. The proportionality of the XLF 
variance to the mean flux holds over more than an order of magnitude of flux 
values. 
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Figure 4: The mean XLF for 14 bright LMXBs from Table 2 (the histogram). 
The confidence level is P(xh > 8-8) « 0.94. 
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Figure 5: The variance of the X-ray flux o as a function of the mean flux 
F x for 14 bright LMXBs from Table 2. The solid line is the best-fit linear 
regression. 
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4 X-ray luminosity finction of LMXB 



The XLF of LMXB is derived from the analysis of observations of point-like 
sources in bulges of galaxies (Gilfanov 2004, Kim and Fabbiano 2004) to have 
the power-law from dN/d\nL x ~ L~ r with 

L x < 2 x 10 37 erg/s 
2 x 10 37 <L x <5x 10 38 erg/s (11) 
L x > 5 x 10 38 erg/s 

Now we wish to show that for systems with orbital periods shorther than ~ 20 
hours the power-law form of XLF can be related to principal evolutionary 
factors driving the mass exchange in LMXBs - the magnetic stellar wind 
from optical star (MSW) and emission of gravitational waves (GW). 




- 1.1 



4.1 Magnetic stellar wind 



Let the mass of the optical star and compact object in a binary be M D and M x , 
respectively, the component mass ratio be q = M Q /M X , and the semi-major 
axis of the binary orbit be a. The characteristic time scale of the orbital 
angular momentum loss due to MSW tmsw = J/ Jmsw is usually calculated 
using the empirical Skumanich law for rotation braking of single late-type 
main-sequence stars, according to which the surface angular velocity of a 
star decreases gradually with age t as \/t. Applying this law to the late-type 
low-mass optical component in a LMXB and assuming the tidal locking of 
the component's spin to the binary orbital period (see e.g. Massevich and 
Tutukov 1988, van den Heuvel 1992 for more detail) yields 

a 5 M x 

™ sw ~ (m x + M D yMt (12) 



In LMXBs the mass transfer is due to the Roche lobe overflow of the optical 
star so we can write 

Ro = af(q) (13) 
where for the shape of the ROche lobe we can adopt 




q < 0.5 
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Then the mean accretion X-ray luminosity is 

M M 5 M x (l + q) 2 f(q) 5 



T~MSW Rl 



Using the mass-radius relation for late-type main-sequence stars R ~ AfJ- 9 ™ 1 
we find 

L ~ M 2 ' 17 - 2 - 67 , a « 2.17.. .2.67 

Substituting this into Eq. (J|]) we obtain for the stationary mass function 
of optical comonents in LMXBs (3 s t ~ 3.52. ..4 assuming the Slapeter intitial 
mass function. Finally, we arrive at the luminosity function in the form: 

' /A L -1.6..-1.13 ? ~ 1-6.. .1.13 (14) 



dliiL, 



o 



4.2 Gravitational radiation 



In the dipole approximation, the time-scale for the angular momentum loss 
in a binary systems reads 

TGR ~ (Mr + M )M M X (15) 
Using the Roche lobe filling condition a = R j f(q) as above we get 

• M MlMl{l + q)f{ q f 

L ~ M ~ ~ — 

Utilizing R ~ M° 9 we find 

L ~ M~ Q,27 ""~ 0,3 , a « -0.27... - 0.3 

Substituting this into Eq. (JHJ) yields /? st ~ 1.05. ..1.08, so finally we obtain 
the luminosity function 

( ' y V' 16 -, r Gfl ~0.16...0.3 (16) 



dlnL 



o 
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4.3 The break of the mean XLF of LMXBs 



The break of the mean XLF of LMXBs is observed at X-ray luminosities 
~ 2 x 10 37 erg/s, which corresponds to the accretion rate on the neutron star 
M ~ lCr 9 M Q /. In our interpretation, the break should correspond to the 
transition from the MSW scale to GW scale as the mass of the optical star 
decreases below a certain value where MSW becomes inefficient. From the 
simle equality tmsw = t gr we find the stellar mass 

M Q « 0.4M© 

which is close to the generally accepted value ~ O.3M (Spruit and Ritter, 
1983). At such masses the mass exchange rate in LMXB drops to 

M(O.4M ) ~3x l(T lo M /yr 

(see, for example, Rappaport et al. 1984 and later calculations of LMXB 
evolution with non-degenerate companions), which roughly agrees with the 
observed value of the break. 



5 Discussion and conclusion 



The above analysis allows us to conclude that the observed mean XLF of 
LMXBs can be generally explained by accretion on neutron stars from Roche- 
filling non-degenerate component driven by gravitational wave emission (be- 
low L x ~ 2 x 10 37 erg/s) and by magnetic stellar wind (above L x > 2 x 10 37 
erg/s). We have shown that these mechanisms, which determines the evo- 
lution of binaries with orbital periods shorther than 15-20 hours, lead to 
power-law dependences of the mass transfer rate M on the optical star mass 
M , which is the necessary condition for the power-law form of the popu- 
lation of such sources. In LMXBs with longer orbital periods, the Roche 
lobe is filling by a subgiant optical star. Using the analytical treatment of 
the mass transfer process in such binaries (Yungelson and Livio 1998), it is 
straightforward to show that in LMXBs with subgiants, too, the power-law 
dependence of the mass exchange rate on the optical star mass is recovered. 

We should also note that a sizable fraction of LMXBs actually recides or 
was born in globular clusters due to effective dynamical interactions between 
stars. The evolution of LMXBs in dense stellar clusters is more complicated 
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than outlined above. So a more accurate evolutionary treatment of XLF 
of such LMXBs is required. For example, the population synthesis method 
seems to be promising. 

The contribution of ultrashort-period binaries dynamically formed in globu- 
lar clusters (like X 1820-30) into XLF of LMXBs was discussed by Bildsten 
and Deloy (20). In these binaries, the Roche lobe is filled by a low-mass 
degenerate white dwarf with inverse mass-radius relation and mass transfer 
is controlled by gravitational wave emission. The analysis of these authors 
shows that a collection of such ultracompact binaries could significantly con- 
tribute to the XLF of LMXBs in the luminosity range 6 x 10 37 < L x < 5 x 10 38 
erg/s and give the slope of XLF as observed. 

Our analysis indicates that the mean XLF of individual sources constructed 
from the RXTE ASM data can be satisfactorily fitted by a quasi-Lorentzian 
distribution with variance proportional to the mean luminosity in a wide 
range of X-ray fluxes. We have shown that such a shape of the mean in- 
dividual XLF does not affect the power-law shape of the XLF formed by 
the entire population of the binary sources. The latter is determined by the 
dependence of the mean X-ray luminosity of each source on the mass of the 
optical star filling the Roche lobe. 

We stress that it is the possibility to describe the mean XLF of individual 
LMXBs by a single law with specific symmetry (proportionality of the vari- 
ance to the mean value) that allows us to recover the power-law form of the 
XLF of the entire source population from evolutionary properties of mass 
transfer in LMXBs. Observations show (Kim and Fabbiano 2004) that XLFs 
of individual galaxies are somewhat different, the unique poer-law of XLF 
is obtained after averaging over many galaxies. This apparently reflects in- 
dividual features of LMXBs in different galaxies (e.g., various contributions 
from ultra-compact binaries, globular clusters, etc.). These observations can 
be used to study formation and evolution of LMXBs in other galaxies. 

The authors acknowledge M.Gilafnov and A.Vikhlinin for useful discussions 
and L.R.Yungelson for comments. The work is partially supported by RFBR 
grants 04-02-16720 and 03-02-16110. The research of PK was supported by 
the Academy of Finland grant 100488. We made use of RXTE ASM data 
publically available from HEASARC of GSFC. 
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